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Abstract 

The chiral phase dependence of fermion partition function in spher- 
ically symmetric U(l) gauge field background is analyzed in two di- 
mensional space-time. A well-defined method to calculated the path 
integral which apply to the continuous fermion spectrum is described. 
The one-to-one correspondence between the nonzero energy continu- 
ous spectra of two pertinent hamiltonians, which are defined by the 
Dirac operator to make the path integral well-defined, is shown to be 
exact. The asymptotic expansion for the chiral phase dependence in 
(m is mass of the fermion.) is proposed and the coefficients in 
the expansion are evaluated up to the next-to-next leading term. Up 
to this order, the chiral phase dependence is given only by the winding 
number of background field and the corrections vanish. 
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1 Introduction 



The quark partition function Z in given gluon background has nontrivial 
dependence on chiral phase a, 

a = arg det m, (1) 

of the mass matrix m |]l], ||. The ^-term in QCD lagrangian comes from 
this a-dependence and leads us to strong CP problem [|l], ^ 

In the previous work we have evaluated the a-dependence, din Z/ da, 
by calculating functional Jacobian in fermionic path integral ||^, ^ . The path 
integral is calculated by solving eigenstates and eigenvalues of two pertinent 
hamiltonians, and H^, which is constructed from the Dirac operator 
of the model. Applying a box quantization to each eigenstates, we discre- 
tised the spectrum and have evaluated the Jacobian in a manifest fashion. 
An important thing which was took for granted in the previous calculation 
was the existence of the one-to-one correspondence between nonzero eigen- 
states for H±: the spectrum of is the same as that of H_ except for the 
zero modes. Under the existence of the correspondence, the contribution to 
din Z/ da from the nonzero modes cancels and only zero mode contribution 
remains. As a result we get only the 6'-term and no other correction terms 
for din Z/ da. 

In a certain background configuration, the box quantization is not nec- 
essarily possible. In this case, we need to start with the infinite size system 
and deal with spectrum with continuous portion from the beginning. The 
one-to-one correspondence then becomes vague. Indeed, there have been var- 
ious intensive studies on this point, especially related to calculation of the 
index of Dirac operators ||, |10|, |11|, [T2|, |T^. It has been pointed out 
that in some models the one-to-one correspondence is no longer hold |Ty, pH]] : 
the spectral densities of H± for the continuous portion are different. If the 
one-to-one correspondence is violated, din Z/ da has extra contribution from 
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the states which belongs to the continuum. 

The purpose of this paper is to develop the calculation of the a-dependence 
in a well-defined fashion which apply for the case that the spectrum contains 
continuous portion and to see if there is any correction in the a-dependence 
of the partition function. We will work in a model of a fermion coupled to 
U(l) gauge field in two dimension. As for chiral property, the fermion 
shares similarity with quarks in QCD: the axial U(l) anomaly is given by the 
winding number density of the background field. Thus we hope the technique 
developed in this paper also apply for QCD. 

We specifically work for spherically-symmetric background field. This 
restriction enables us to analyze the problem in terms of the scattering theory 
for a particle in a centrifugal potential in two dimension. The spectrum is 
continuous and extends to zero energy. We consider the background as having 
any real winding number 



which prevents us from compactifying the space-time and from discretising 
the spectrum. Depending on whether the winding number is zero or not, 
the potential possesses 1/r^-tail at spatial infinity. This fact causes non- 
analyticity of the scattering waves at zero energy and will turn out essential 
for nontrivial din Z/ da. 

This paper is organized as follows. In section 2, we review the calcula- 
tion for din Z/ da for discrete spectrum and rewrite it in terms of spectral 
density. This form of dlnZ/a enables us to calculate it for the continuous 
spectrum. In section 3, we describe the definitions of spectral densities and 
write din Z/ da in two terms, a contribution from nonzero energy continu- 
ous spectrum and those from zero energy states. We use contour integral to 
regularize the zero energy contribution. In section 4, we define Jost func- 
tions in two dimensional scattering problem and write the spectral density 
in terms of them. In section 5, we show the one-to-one correspondence for 




(2) 
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the nonzero energy states is exact. In section 6, we will evaluate the zero 
energy contributions. We adopt the asymptotic expansion of din Z/ da in the 
inverse of squared mass parameter l/|mp and calculate the coefficients up 
to the second order term. Section 7 is devoted for summary and discussion. 

2 The a-dependence of Z 

We take a fermion ip in two dimension with a complex mass m and coupling 
to external U(l) gauge field A^. The lagrangian in Euclidean metric is of the 
form 



(9^ - ieAu_)ai^ - m 



1 + 



— m 



I -a, 



(3) 



2 7 V 2 

where e is the coupling constant, (/i = 1,2) are the first and second 
components of Pauli matrices. 



1 

the third component is used to define chirality. 

The chiral phase a is simply related to m as 

m = Imle*". 



-i 

1 



1 
-1 



(4) 



(5) 



(6) 



We intend to examine the a-dependence of the partition function 

Z = j[dilj][di)]e-S'^''. (7) 

Specifically we evaluate the derivative din. Z/ da. By definition, it is given by 

dlnZ 1 fr , r , T, /■ ,o / dL\ 

da I 



da 



= ^ /[#][#] / 



d\\ -''-^\ e- ! 



= j d\ l^^ 



m 



1 + CT, 



m 



(8) 
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where (O) denotes the expectation value of operator O in a given gauge field 
background. 

The evaluation of Z and dXwZjda in terms of fermionic eigenstates in 
gauge field background has been given by discretising the spectrum with box 
quantization P] . Let us briefiy review the results. The Dirac operator L has 
an explicit form 



with 



D = i{di-ieAi) + {d2-ieA2), (10) 
= ^{^l-^eA^)-{^2-^eA2), (11) 

for lagrangian (^ . It is not hermitian for complex m and thus we cannot use 
eigen modes of C to define the functional space. Instead, we use eigenstates 
0^'' and (j)^^'' of "hamiltonians" 

= DD^ = -{d^-teA^f-eFu, (12) 
H_ = D^D = -{d^-teA^f + eF^2, (13) 

where n is integer and used to label different eigenstates. The eigenvalues 
are non-negative. The spectrum of H± coincide for nonzero eigenvalues. 

This can be seen by noticing the one-to-one correspondence between and 

An) 

r- , 



D^<P':\ (14) 



0?) ^ -±=D0, (15) 

where E^") = ^ = E^^\ 

Expanding the field ip and ip in terms of these eigenstates, diagonalis- 
ing the action by an variable transformation for Grassmann path-integral 
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variable, and taking into account the corresponding Jacobian, we get 



J\m 



n_|_+n_ 



(n) 



+ \vrt\ 



(16) 



where n± are the number of zero modes for li± and the prime on the prod- 
uct means the product must be taken over each pair of nonzero eigenvalue. 
J denotes the Jacobian in the path-integral measure and it carries the a- 
dependence as e*"("+~"-). Further the expectation value in (H) is given by 



cilnZ 

da 



E 



\m\ 



(n) 



(17) 



Because of the one-to-one correspondence (p!^)-([T5|), the nonzero eigenvalues 
cancel out and (jl^) is consistent with the a-dependence of J . 

We extend these results to the case of continuous spectrum by generalizing 
them in terms of the spectral density, the number density of states per energy. 
For the discretised spectrum, it is obviously given by 



We rewrite (|T^) and ([17|) as 

roo 1 

lnZ = lnJ^+ / rfE-[p+(^) + p_(^)]ln(E + |mn 



and 



cilnZ 

da 



dE 





m 


2 


E + 


m 


2 



[p+{E) - p^{E)] 



(18) 

(19) 
(20) 



We take (|1^) and ( pOD as definitions for Z and its a-dependence in terms of 
spectral density p±] we develop the way to calculate p± that is applicable to 
the continuous spectrum. 

We need to mention two problems about the definitions (p!^)-pO|). First, 
we have to be careful in the ii^-integrals at ii^ = to get well-defined results 
when the continuous portion in the spectrum extends to zero energy. A more 
careful definition which prevent the problem will be given in the next section. 
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Secondly, the spectral density behaves as p± (E) ~ V^/47r as — > oo with 
the volume V of quantization box in two dimension. Thus the E'-integral 
in In Z has UV divergence. We can define once-subtracted spectral densities 
6p±{E) which is defined as the deviation of p{E) from the spectral density 
V/An of the free hamiltonian. In InZ we use 6p±{E) instead of p±{E). The 
perturbative calculation for InZ, which measures its deviation from the case 
of free fermion, reveals that the correction to InZ is finite [|l^. Thus the 
one subtraction suffices to make p±{E) and \nZ finite. As for din Z/ da, the 
divergence cancels between p+ and p_ and no ill-definedness takes place. 



3 Spectral density 



The problem is now to develop a decent method for evaluation of the spec- 
tral densities p±{E). The various methods have been intensively studied in 
connection with fermion number fractionization and with Witten index in 



supersymmetric quantum mechanics [11, 12, 13, 14 . A subtlety of the cal- 



culations related to the order of various integrations are emphasized for the 



Dirac operators whose continuous portion extend to zero energy [|TT|. Here, 
we adopt a different way of defining spectral densities. The advantage of our 
description is that it naturally leads us how to regularize the contribution 
from the zero energy states. 

We start with heat kernel K{t; x, y) for H. (We suppress the subscripts 
± until we need to distinguish H±.) It is defined by the equation 



K{t- x, = for r > (21) 
and the initial condition 

X(0;x,y)=5(x-y), (22) 

where 5{x — y) denotes the two dimensional Dirac (5-function. A set of (pl|) 
and (E2) fixes K{t; x, y) completely. 
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The resolvent R{E; x, y) is then defined as the Laplace transform of the 
heat kernel K{t; x, y) by 

roo 

R{E-x,y)= / dre^^K^T-.x^y). (23) 

Using the definitions (pl|)-(|23|), we see R{E;x,y) satisfies 

{H~E)R{E-x,y)=6{x-y) (24) 

as long as the integral in (|23| ) converges. According to general properties of 
Laplace transformation, R{E; x, y) for fixed x and y is a. holomorphic function 
of complex E whose real part is less than some critical value Ec, (Re E < Ec). 
In our case, the spectrum of H is nonnegative thus Ec = 0. We analytically 
continue R{E; x, y) in a complex E-plane with a cut along the real positive 
E-ax.is. 

Obviously, K{T;x,y) is restored by the inverse Laplace transformation 
K{t; x,y) = ^^ dEe-^'-RiE; x, y), (25) 

where the integral contour C runs from c — too to c + ioo, (c < 0). 

We define a state density p{E; x) of H, the number of states per unit 
energy and unit space volume, by the discontinuity of the resolvent R{E; x, y) 
at the cut |jT6|, 

piE] x) = lim lim — \R(E + ie; x, y) - R(E - ie; x, y)] . (26) 

The order of two limitations in (^), (e +0) before {y —>■ x), cannot be 
altered. This circumvents ill-definedness which happens in the space with 
dimensionality larger than one. In appendix ^, we give simple examples for 
free hamiltonians in various dimensions which show the above definition is 
consistent with our understandings to p{E] x). 

The spectral density p{E) is then defined formally by 

p{E) = J d\p{E;x), (27) 
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and the once-subtracted spectral density by 

5p{E)^ j d\\p{E-x)-^, (28) 

where is the state density for the free hamiltonian. This subtraction 

make the ii^-integral in ([T^) finite. 

Now let us turn into a regularization problem at = 0. In the background 
we consider, the resolvents are singular there. There are two reasons for the 
singularity: the existence of normalizable zero modes 0'-°^ adds pole-type 
singularity 

i?(i?;a;,y)~i0(°)(x)0(°)t(^) (39) 

to the resolvents; since the continuous portion in the spectrum extends to 
zero energy, the resolvents may have other type of singularities as well. The 
definition ( p6D is not sufficient to take into account these singularities. 

To answer this problem, let us recall the calculation of the anomaly, 
the anomalous term in Ward-Takahashi identity. According to Fujikawa, it 
is given by taking a functional trace of chirality cXc using the heat kernel 
regularization 0. It corresponds to r — > limit of 

A{t,x) = lim [K+{T;x,y) - K_{T;x,y)] . (30) 

y >x 

[The limit [y — ^ x) should be taken before (r —>■ +0), otherwise we get 
trivially zero as the result because of the initial condition (^).] This quantity 
lim^_+o -^(t, x) must be the difference between the state densities per unit 
spatial volume for Hj^ and we use eigenstates of H± for the basis in 
the functional space and they have definite chirality ±; thus the difference 
between the state densities generates nonzero Jacobian in the path integral 
measure under a chiral transformation e^^^^^'^" 0]. We use the Laplace 



transform (25) for the heat kernels in ( pOD and bend the integral contour C 



so that it pinches real positive E axis leaving a small circle C around the 
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origin, to get 

AiT,x) = \im—idEe-^^[R+{E;x,y)-R.{E-x,y)] 

/oo 
c/Ee-^"[p+(E;x) -p_(E;x)], (31) 

where C is a clockwise oriented with radius e. (See Fig. 1.) The second 
term represents what it must represent: it is the difference between the state 
densities contributed from the nonzero energy states. The first term tells us 
how to regularize the contributions from zero energy states. If the singularity 
comes from normalizable zero modes, the contour integral over C properly 
gives the density by zero modes. The expression in (|3T|) , however, applies to 
wider class of singularities. 

Incidentally, note that the calculation of the anomaly in term of A{t, x) 
is gauge invariant procedure: two different heat kernels for and its gauge 
transform A^ = A^ — (9^A / e are related by 

K{t; X, y) = e^(^(-)-^(2^))ir'(r; x, y); (32) 

the difference disappears after we take the limit (?/ — > x). 

Applying the same reason to regularize the contributions from zero energy 
states, we finally adopt the definition^ 

d In Z c X y I ttIj | ^ 

+1 r ^^^TT-^ t-^^^^) - '^-(^)] • (33) 

4 Jost function in two dimension 

We now apply the method described for spectral density in the previous sec- 
tion to the case of spherically symmetric background field. The problem 
similarly regulated expression for h\Z in terms of Sp±{E) is easily given. 
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reduces to solving Shrodinger equations of two dimensional scattering. Ac- 



cording to a well-established technique in three dimensional scattering [|T| 
we define two types of solutions for the Schrodinger equations, regular solu- 
tions and Jost solutions, and define Jost function in terms of them. We can 
derive resolvents and spectral densities in terms of them. Especially we will 
see that spectral density is related to the Jost function. 
Let us define the polar coordinates r and by 

a;i = rcos6', a;2 = 'rsin6'. (34) 

We take radial gauge Ar = Q and choose Ae which depends only on r, where 

Ar = Aicos6 + A2sme, (35) 
Ae = -AirsinO + A2rcos9. (36) 

The resulting field strength = dAe/rdr is then spherically symmetric. 

The Ag must vanish at the origin so that is continuous. To avoid any 
subtle singularity which may appear when we calculate derivatives of with 
respect to coordinates at r = 0, we further restrict Ag to such a configuration 
that its derivatives of any order with respect to r vanish at r = 0. An simple 
example of such a configuration is Ag ~ e~^^^. 

For the background to have nonzero winding number, Ag approaches to 
a constant, 

eAg ~ z/, (37) 

as r goes to infinity. Note the constant u is the winding number in (|^). 
In this background, the derivative operators (|TU|)-(|n|) are 

i^dr r 06 r ) ' ^^^^ 
^'•--•"f^+Z^+^V (39) 
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The Schrodinger equations H±(f)± = E(f)± with E = k"^ read 
(P 1 



W 2 + ^(^'-7)+^±W 



(40) 



where '?/'"(r) and '?/'"(r) denote partial waves with definite integer angular 
momentum n 



1 



and the potentials are 



V+(r) 



e^Al - 2neAo] - (± 



1 d{eAe 



r dr 



(41) 



(42) 



Since the procedure hereafter applies equally to both and ip'^ , we suppress 
the subscript ±. 

In background free case, V{r) = 0, the typical two solutions of (|40| ) are 



Jn{kr) for n > 



;-)'"7y'^h(^^) forn<0 



and 



w'^{k,r) 



ixlYH^^K^^r) forn>0 

)l"^/^ifW(A;r) forn<0 



(43) 




(44) 



where J„ and H^^ denote the Bessel function and the first kind Hankel 
function of order n, respectively. The phase convention for n < in (^3])- 
( ^if ) helps expressions below to be simpler. (See, e.g., (|H^)-(p3|).) The free 
solutions asymptotically approximate to 



u"'(k,r) ~ —= cos \kr 
Vk \ 



{n + l/2y 



w'^{k, r) 



Vk 



exp 



i \ kr — 



(45) 
(46) 
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and their Wronskian is 

W [w'^ik, r), vJ'ik, r)} = w'^ik, r)^^^-^tlll _ ''^ ^"(A;, r) = 1. (47) 

dr dr 

Based on these free solutions, we define regular solution ip'^{k,r) and Jost 
solution f"'{k,r) by integral equations 



r 



i^'^ik.r) = u'^ik.r)- dr'\w''{k,r)u''{k,r') 

Jo 

- w"(A;, r')M"(A;, r)] F(r')<^"(A;, r') (48) 



and 

r(fc,r) = w"(A;,r)+ / dr' [w''ik,r)u''{k,r') 

Jr 

- w^'ik, r')u''{k, r)] V{r')r{k, r'). (49) 
These equations are solved by iteration and defines a series for both and 

The advantage of defining the solutions by the integral equations, or 
the series expansions generated by the equations, is that we can regard the 
solutions as analytic function of complex values of k for fixed real positive r. 
By the same estimation about convergence of the series as in Ref. we 



can see the following analytic nature of the solutions. For (p"'{k, r), the series 
absolutely converges if the potential meets a condition 

• r'\V(r')\ 

dr' ' , < oo. (50) 

l + \k\r' ^ ' 

Then, since each term in the series is an entire function of fc, <^"(A;,r) is 
also an entire function of k. [Strictly speaking the condition (|50| ) is for 
n 7^ 0. Note that H^\kr) includes a term proportional to Jn{kr) Inkr. For 
n = 0, this term requires that the potential obeys a slightly stronger condition 
Jq dr'r' \V{r')\ Inr' < oo at r ~ 0. Note also that the part proportional to 
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In k cancel in (^) and does not affect the entireness of tlie series with respect 
to k.] Similarly, the series for ^^{k^r) absolutely converges if ImA; > and 

■*'™<oo. (51) 

f"'{k,r) is, then, a holomorphic function in upper half /c-plane. 

Note the potential ( ^21) with (^) has 1/r^-tail at infinity for nonzero 
winding number u. Thus the point /c = is the singular point for the Jost 
solutions. 

The regular solution and Jost solution satisfy some symmetry properties. 
They are derived by noticing that [<f^{k, r)]* and (f^{—k*,r) satisfy the same 
wave equation and are both regular at r ~ 0. Thus they should be propor- 
tional to each other. The same reason applies to ip"'{k,r) and ip"'{—k,r). By 
using the behavior at r ~ 0, 

^"(A;,r)~M"(fc,r)~^fcHy , (52) 

we obtain 

[<^"(fc,r)]* = (-)l'^l(^"(-r,r) and (^"(/t, r) = r). (53) 

Similarly, noticing /"(A;,r) ~ w"{k,r) at r ~ 00 and (ffBI), we get 



[r{k,r)r = {-pr{-k*,r). (54) 

We define the Jost function in terms of Wronskian by 

F"(fc) = iy{r(fc,r),^"(A;,r)}. (55) 
F"{k) is a holomorphic function in upper fc-plane and obeys a condition 

F"(fc) ^1 as |A;| ^ oo (56) 
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for the potentials that meets condition ( |50D with r = oo [|T^. Further by 
(H) and fl), we see 

= (57) 

ip^{k,r) can be expressed in a linear combination of f^{k,r) and a solution 
that is linearly independent from f^{k,r), say g'^{k^r). The Jost function 
represents an amplitude of ip"'{k,r) with respect to g"'{k,r) in the linear 
combination. Especially, for real positive k, f'^{k,r) and ^^{—k^r) are the 
linearly independent two solutions. The linear combination for (f"-{k,r) in 
this case is given by 



<^"(fc, r) = - [(-)I"IF"(A;)/"(-A;, r) - F''{-k)r{k, r)\ (real positive A;). 

^ (58) 
Now we can write the resolvent R{E;x,y) using (f"'{k,r), /"(A;,r), and 
F"'{k) in the following way. First define a function R^{k\ r, r') which is holo- 
morphic in the upper half /c-plane by 

1 



R"{k;r,r') ^ ——[e{r-r')r{k,r)y,"{k,r') 
+e{r' -r)<^''{k,r)r{k,r')]. 



(59) 



R"'{k;r,r') satisfies 



n--]+ V{r-) - r 



i?"(A:;r,r') = 5(r-r'). (60) 



Then R{E\ x, y) for E = k"^ is given by 



R{E;x,y) 



5 '"yjj 



(61) 



where {vy) is the polar coordinates for x (y). The resolvent defined this 
way obeys ( p3D and is holomorphic in complex ii^-plane with a cut along the 
positive real axis. 
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One may wonder an ambiguity to add an linear combination of terms 
made of the product of the solutions, like ip^{k,r)f"'{k,r'), to R^{k;r,r'). 
It does not change the equation (|60D. Two requirements fix this ambiguity: 
R{E; X, y) must be regular at both r^. ~ and ^-^^ ~ since the background 

is regular there; R{E;x,y) must be a decreasing function of \x — y\ as 
\x — y\ CO so is K{t; x, y). 

Remembering the definition for the state density (^) and using the rela- 
tion (|58D for real positive fc, we get 



1 



2 



P(E;x) = —- E ' ^.r... (62) 



We further calculate the once-subtracted spectral density 6p{E). In ap- 
pendix 0, we show Sp{E) is related to the Jost function hyf\ 



-oo 



5 Symmetry for Jost functions 

So far we have defined regular solutions ip^ and Jost solutions /j! and ex- 
pressed the spectral density by Jost functions. In this section, we investigate 
about the one-to-one correspondence for these solutions and show the con- 
tribution to din Z/ da from the nonzero states vanish. 

We first prove a set of relations which hold for the solutions; 



^The relation can be given by using the equivalence of Jost function and Fredholm 
determinant for the case of H± being Fredholm. In the problem we are considering, H± 
are not Fredholm operators. 
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where ip^ denotes both and (p^ collectively. These relations are the partial 
wave versions of the relations (|T^)-(|T5[); note the differential operators D 
and D^, (]38|)-(|39[), connect only the partial waves whose angular momentum 
differs by one. Also, these relations holds for any complex values of k at which 
the solutions are finite and well-defined. In this sense, they are extension of 



(p!4D-(|T5|) to complex values of energy. 

We specifically prove (|6^ for (p^{k,r) and y9"^"^(fc,r); the others can be 

shown in quite similar way. Let us use notations 

^+ d n + 1/2 eAo ^+ eAo 

= — - + - = T>1 (66) 

dr r r r 

_ d n + 1/2 eAg ^ eAg 

V = — + -=Vo -. 67) 

dr r r r 

Note a set of formula, 

yVlu''(k,r) = u''+\k,r), lVou''+\k,r) =u''(k,r), 

k k 

yVlw''(k,r) = w"+Vfc,r), ^VQw''^\k,r) =w''(k,r), (68) 

k k 

which apply to the free solutions defined by (|I3|)-(P^. The equations (^) 
for ip^[k,r) and ip^~^^{k,r) in abbreviated but manifest notations are 

k^ - VqVI) 
k^ - VIPq 

The potentials V+ with the angular momentum n and with n + 1 are also 
written as 

V+ = -VqA-AvI + A^ (71) 
V_ = -VlA- AVo + A\ (72) 

where A = eAg/r. Using explicit integral expression for 1/(A;2 - VoVl) m 
terms of and and applying the formulas (|68D , we verify 



Vl r = ; — Vl (73) 
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as long as it operates on a function, say h{r), that is regular enough at r = 
so that 

limw''+\k,r)h{r) =0. (74) 



This condition is satisfied for the potentials derived by ( ^21 ) with Aq that 
satisfies the regularity condition at r = which we have mentioned in section 
4. Similarly, we get 



1 



Vn. 



Now applying Vq to (|69D and using ([73|) , we get 



where we have used 



Vlu 



vM+[k''-VlVo)A 

2 



A;2 - vlV^ 
A {k^ - VoVl) - A^Vl + VlA 



Thus obeys 



1 



A^dJ - VIA 



t 42 



fc2 - vlV, 

Iterate 93+ in the second term by (|69| ) and get 

1 



+ 



A:2 - ptp, 



-V^ku 



0^0 



A;2 - VlVo k^ - VlVo 



k^ - VlVo k^ - VlVo 



V.Vl-A [k^ - VoV, 

A'vi-vm 



(75) 



(76) 



(77) 



(78) 



(79) 
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Continue the similar iterations and get 

oo / 



m=i \k^-VlVo 



V- ku 



n+l 



(80) 



This series coincides what is generated by (^) up to a factor fc; thus we get 

The proof for f^{k,r) and f!^~^^{k,r) is quite similar. The equations (^) 
and (|7^) holds as long as they operate on a function h{r) that obeys 



lim M"+^(A;,r)/i(r) = 0, 



(81) 



which is also satisfied in the problem we are considering. 

The relation (|6^)-(|65D lead very critical result. The contribution of the 
nonzero energy states to din Z/ da is given by the difference between the 
spectral densities of if+ and According to (|63|) , it reads 



with k 



F^+\k) FI^+\-k)_ 
E. On the other hand, because of (|64D-(|f^ 



we see 



W{f'l+\k,r),cp-_-'\k,r)} 

1 [v^flik,r)VV^^lik,r)~VV^flik,r)V^cpl{k,r] 



(82) 



W{fl{k,r),^lik,r)} 



I.e. 



F!^+\k) = Fl{k). 



(83) 
(84) 



Thus the contribution of nonzero energy states to d\n.Z/da vanish. 



6 Contribution around E = {) 

We have seen that d\nZ/da gets contribution only from the zero energy 
states, i.e., the integral over C, (see Eq. (pBf)). To know how the resolvents 
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behaves at ~ is essential to evaluate the integral. Due to the 
tail in the potential, the analysis is quite difficult. Instead, we attempt to 
estimate din Z/ da by an expansion in l/|mp: the i?-integral is now carried 
out over finite region; thus we expect the expansion of \m\'^/{E+ |mp) in the 
integrand in l/|mp results in a decent expansion, which is at least asymptotic 
to d In Z/ da as |m| goes to infinity. Let us define 

Sk= /Alim-^i" dE{-E)''[R+{E;x,y) - R_{E;x,y)] (85) 

and assume they are finite. The expansion is then given by 
d\nZ ./ 1 ^ 1 



i\So + --,Si + -—S^ + ...\. (86) 
da \ \m\'^ |mp j 

In order to get the coefficients Sk-, we look into 

A(r) = / d\A{T,x). (87) 



Recalling Eq. (0), we notice that A(r) also get the contribution only from 
integral over C and thus has the expansion in r with the coefficients Sk- The 
advantage of considering A(r) is that the asymptotic expansion of ^(r, x), 

A{t, x) ~ Aq{x) + tAi{x) + t'^A2{x) + (88) 

can be evaluated by a straight forward way^; once we get Ak{x), Sk are given 
by integrating them, 

Sk = k\ f d^xAkix). (89) 



We evaluate the coefficients Ak{x) using the explicit expression for A{t, x) 
Ar, X) = Jim / ^ (e--^- - e"^^") e-(^-). (90) 



■'The interplay between the index of Fredholm operator and the asymptotic expansion 
of its heat kernel is described in 
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We regard the operation of e""^^* is calculated by the corresponding power 
series. To take the limit y —>■ x, we move the factor e*^'^ on the right side 
of the operators to left. Every time one commutes the partial derivative 
with e*P'^, we changes it with + ip^. Scaling the momentum by i/r, we 
get 



A{t, x) 



1 r 2 ^ ]_ 

r J {27ry n\ 



T 



2ip ■ id — ieA) ,^ „ 



-(Fi2 ^ (91) 

The expansion (|88|) is given by picking up terms with the same powers in r. 
We have evaluated the series explicitly up to A2 and got 

Aq{x) = 7^eFi2, 
Ztt 



1 



A2ix) = ^dldleF,,. (92) 

Note all these terms must be gauge invariant as well as having parity odd. 
It is further interesting to note that the results are total divergence: even 
under the requirement of the gauge invariance and the parity, the term pro- 
portional to (-^12)^ was possible for ^3; these terms cancel out in the explicit 
calculation. 

The coefficients Sk are now calculated by (|55|). The leading term is 



5o = (93) 

i.e., the winding number of background field. This is exactly the same way 
as in four dimensional QCD. Since F12 approaches zero faster than 1/r^, the 
second and third correction terms are zero. 
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7 Summary and Discussion 



In summarizing, we have evaluated the chiral phase dependence of free energy 
In Z in two dimensional spherically symmetric background. We divided up 
the contribution into two; one from the nonzero energy states and one from 
zero energy states regularizing the later by the contour integral over C. 
The contribution from the nonzero energy states is shown to be zero. In 
other words, the one-to-one correspondence is exact. As for the zero energy 
states contribution, we proposed the expansion in l/|mp. The coefficients 
are given explicitly by the behavior of the resolvents at ~ 0, see Eq. (|85|). 
We have evaluated the coefficients up to the second power of l/|mp. We 
found the a-dependence of the free energy is given by the winding number 
of the background U(l) gauge field and corrections to this relation is zero up 
to l/|m|^. 

The expansion for Air^ x) also implies that the expectation value of the 
axial U(l) current can be given by the expansion of l/|mp. Comparing our 
result with the anomaly equation for the axial U(l) current 

it 

TT 

we get 

67r|mp 307r|m|^ 



— — ^6 / — \ 

m'4){l + (Tc)ip - m*ip{l - (Tc)ip = —Fi2 - df, itlja^acip) , (94) 

TT ^ ' 



Ua,a^^) ~ (-.) + + ...). (95) 



Finally, it is worth mentioning the interplay between nonzero winding 
number and 1/r^-tail in the potential. The existence of the 1/r^-tail means 
the violation of the Levinson theorem in the scattering problem [0. The 
violation indicates that the total number of states in the spectra of H± are not 
the same. This difference causes the nontrivial a-dependence of the Jacobian 
in the path integral measure when we rotate a away from the mass term in 
the lagrangian. 

Acknowledgement 



22 



This work is supported in part by the U.S. Department of Energy under 
contract No. DE-AT03-87ER40327. 



A Appendix 



We verify in this appendix that the definition ( P6|) for the state density 
p{E; x) gives correct resuhs for hamiltonians which describe a particle moving 
freely in d-dimensional space. The hamiltonian is 

V 

dx,. 



H = -(^] (96) 



with an unit 2m = 1 for the particle mass m; fi runs from 1 to d. According 
to the standard box quantization, the number of quantum states per unit 
momenta and coordinates, is 

= ^ (97) 

U^dp^dx^ (271)^ ^ ^ 

in an unit h = 1. The state density p{E]x), the number of states per unit 
energy and unit coordinates is thus 

where s{d) is volume of unit sphere S'^^^, s(l) = 2, s(2) = 27r, s(3) = 47r, ... 
(s(l) is the number of points that satisfy p"^ = E.) 
For d = 1, the resolvent R{E; x, y) is 

R{E-x,y) = ^^^e-^\^-y\ (99) 

where we choose a cut of \/—E along positive real i?-axis. Then, 

piE- x) = lim — ( , ^ , ^ I = (100) 
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Similarly, for d = 2, 



1 



R{E;x,y) —\n(V^\x-y\) for |s - y| < (101) 



and 



p{E-x) 



For d = 3, 



and 



lim 

e^+o 2iTi 

1 

An' 



-ln(v/=E37i)+_ln(v/=E + 



27r 



le 



R{E-x,y) 



4:7r\x — y\ 



-E\x-y\ 



(102) 
(103) 



p{E; x) = lim lim 



y^x 2m AtcIx — y\ 
E 



^-V-E-ie\x-y\ _ ^-V-E+it\x-y\ 



47]-2 



(104) 



All these results are consistent with (p8|). Note that for d > 2, R{E;x,y) 
diverges as y ^ x. Thus in order to avoid the ill-definedness, we have to 
calculate discontinuities at the cut first before taking the limit y ^ x. 



B Appendix 

In this appendix we show (|63|) . To see this, we calculate 



dk 



+ 



dk 



F"(A;) F'^{-k*) 
F^ + F^{-k*) 



and take a limit that k approaches to real positive values, k ~ |A;| +ie, where 
'dot' denotes derivative with respect to k. Using the definition of the Jost 
function (ESl), we obtain 



I{k) = h{k,r) + h{k,r), 



(106) 
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where 



and 



h{k,r) ^ -l-^W{r{k,r),^-{k,r)} 

+ Fr^{-k*) ^ {r(-^*, r), r)} (107) 



h{k,r) ^ -l-^w{r{k,r),^-{k,r)} 



F'^i-k 



]_.W {r(-F, r), ^"(-r , r)} . (108) 



Since 



0-(_r,r) = (-)H+i[0"(A;,r)r (109) 
and ip"'{k,r) becomes real as k approaches to the real value, 

2i 

h{k,r) ^ {y-(fc,r),0"(fc,r)} as 6 ^ 0, (110) 

where we have used (^) and (|58D. Noting that ip"-{k,r) obeys an equation 
^2 1 / 2 r 



+ ^ - - + F(r 



0"(A;, r) = A;V"(A;, r) + 2A;(^'^(fc, r), (111) 



which is given by differentiating (|40|), we obtain 

{¥."(A;, r), vj"(A;, r)} = -2A; [y."(fc, r)]^ . (112) 
Integrating this equation with the boundary condition 

limiy {(^"(A;,r),0"(A;,r)} = (113) 



(see equation (|52|) ), we get 



/i(A:,r) -.—————/ rfr'[^'^(fc,rO] as e ^ 0. (114) 
^ ' ^ FHk)F'' -k) Jo ^ ^ 



F^{k)F'^{-k) Jo 
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On the other hand, l2{k,r) becomes independent from specific potential 
of the scattering problem in the limit of r — > cxd, which can be seen in the 
following way. /"(A;,r) approaches to the free solution and thus 

r{k, r) ^ w'^ik, r) as r ^ oo. (115) 

ip"'{k,r) consists of exponentially decreasing component e~"' and increasing 
component e'^'" at the complex value of k. Since f'^{k^r) decreases exponen- 
tially, only the increasing component can remain in ( |108| ) as r ^ oo. The 
Jost function is the coefficient of the increasing component in the expansion 
of (p"'{k,r). Because of (0), we can specifically write it as 

r) F'^{k)u^{k, r) as r oo, (116) 

where we are only interested in the increasing component in u"'{k,r). Plug- 
ging ( |115D and ( |116|) into (|108|) , we see that the potential dependence, the 
Jost function F"'{k), cancels in numerator and denominator. 

I{k) for real k is now written by taking r — > oo in Ii{k, r) and l2{k, r) as 

^C^) = }Sk { F-(k)tH-k) [ *■' '-'M' + '■') • 

and the contribution of the second term is potential-independent. Note the 
first term is related to the integral of p{E; x) over x. Further note F'^{k) = 
and I{k) = for the free potential; thus /2(/c,r) should satisfy 

hik.r) ^ -Aik r dr'[u''{k,r')f as r ^ oo, (118) 
Jo 

which is related to spectral density for free potential. Thus Sp{E), the de- 
viation of spectral density from one for free potential, is given by summing 
I{k) over all angular momentum n as 

-| oo 

SpiE) = ^ E Hk), (119) 

which is (|63D. 
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Figure captions 

Fig. 1. The integral contour C on the complex plane with a cut along 
the real positive axis. 
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